function [u,l,xx]=gausstype2(x,P)

sigma=P(:,1);
m1=min([P(:,2),P(:,3)]);
m2=max([P(:,3),P(:,2)]);

L=length(x);

u=[];
l=[];

parfor i=1:L
if x(i)>=m1 & x(i)<=(m1+m2)/2
mu1=1; 
% mu2=gaussmf(x(i),[sigma,m2]);
mu2=exp(-(x(i)-m2).^2*(sigma).^2);
elseif  x(i)>(m1+m2)/2 & x(i)<=m2
mu1=1;
% mu2=gaussmf(x(i),[sigma,m1]);
mu2=exp(-(x(i)-m1).^2*(sigma).^2);
else
% mu1=gaussmf(x(i),[sigma,m1]);
mu1=exp(-(x(i)-m1).^2*(sigma).^2);
% mu2=gaussmf(x(i),[sigma,m2]);
mu2=exp(-(x(i)-m2).^2*(sigma).^2);
end

u=[u,max([mu1,mu2])];
l=[l,min([mu1,mu2])];
end

xx=x;
